syms a0 a1 a2 a3
eye(2);
s1=[0 1; 1 0];
s2=[0 -i; i 0];
s3=[1 0; 0 -1];
ss=a0*eye(2)+a1*s1+a2*s2+a3*s3
